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\/  INTRODUCTION 

*  r 

This  paper  concerns  the  application  of  lattice  filtering^3r^73] 

techniques  to  discrete  frequency  domain  data.  Time  domain  models  for  two 

different  antennas  are  found,  the  monopole  and  cavity-back  spiral.  Lattice 

cl 

models  and  a  layered  medium  model^£4,5]  are  obtained. 

The  following  is  a  brief  list  of  important  notation  used  in  the  work, 
basically  that  of  Parker  [1], 


Number  time  domain  sampled  points 
Discrete  time  domain  signal 
Forward  prediction  error 
Backward  prediction  error 
Time  domain  sampling  period 
Discrete  Fourier  transform  of  y^ 
Autocorrelation  coefficient  of  lag  i  T$ 
Reflection  (partial  correlation)  coefficient 
Fundamental  angular  frequency 


Order  of  difference  equation  for  y^  s 
ith  coefficient  in  difference  equation  for  y^  s 
Right  going  planar  wave 
Left  going  planar  wave 
Wave  propagation  speed  in  boundary  1 
Characteristic  impedance  in  boundary  1 


□  □ 


PRODUCTION  OF  AUTOCORRELATION  COEFFICIENTS 

We  assume  all  i«0,...,m  are  known.  The  method  by  which  these 

complex  quantities  are  obtained  is  immaterial.  However,  for  this  work  they 
were  obtained  through  experiments  using  the  EMP  Engineering  Research  Omni¬ 
directional  Radiator  (EMPEROR),  the  heart  of  which  is  an  HP8510  Network 
Analyzer. 

The  calculation  of  R^  i«0,...,m  is  straight  forward. 

Ri  -  | YQ | 2  +  |  l  |Ya|2  cost  <i)gi  Ts  i-0, . . . ,m.  (1) 

ojg  is  related  to  the  sampling  period  in  equation  (1)  by 

“o  *  pnrfi^Ts  (2) 

where  [2m  +  l]2p  is  the  nearest  power  2  operator  on  2m+l.  The  number  (N)  of 
time  samples  is  taken  to  be  £2m+l]2p.  is  useful  in  determining  an  estimate 
of  y.  i-0,...,N-l. 

PRODUCTION  OF  LATTICE  STRUCTURES 

For  a  linear  causal  stationary  real  autoregressive  (AR)  process  of  order 
m,  the  forward  prediction  error  is  given  by 


m 


mn 


"  yn  7^.  ^mt  ^n-l 
*■1 


n  -  0 . N-l.  (3) 
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i*  1  * 


The  backward  prediction  error  (sometimes  called  postdiction  error)  is  the 
following  expression: 

m 

emn  "  yn-m  bml  yn-m+l  n  "  °»**,»  N_1* 

The  bm^s  in  equations  (3)  and  (4)  are  the  difference  equation  coefficient  for 

an  AR  process  of  order  m.  It  is  assumed  that  em„  and  £  have  the  same 
r  mn  mn 

statistical  properties,  that  5„  -  e^  -  y  .  and  that  m  i  250. 

r  f  *  on  on  ■'n' 

Using  the  Levinson-Durbin  algorithm  and  equations  (3)  and  (4)  we  can 
write  the  following  important  relations: 


em+l,n  "  emn  ”  ^m+l  em,n-l 


em+l,n  "  “km+l  emn  +  em,n-l 


* 


where  the  quantity  km+^  (reflection  coefficient)  is  given  by  either  the 


expectation 


E ( e  5m  n  i) 

1  mn  m , n- l 

E<IVn-ilZ) 


ratio  or  equivalently 


km+l  "  bm+l,m+l  * 


Another  way  of  calculating  the  k-values  is  to  use  the  Shur  algorithm  [6], 
since  by  now  the  R^s  are  known.  This  algorithm  Involves  the  following  steps: 
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1.  Define  a  generator  matrix 


R0  R1  •**  Rm 


0  R<  •  •  •  R_ 
1  m 


2.  Shift  the  first  column  of  the  generator  matrix  and  define 


0  R0  Rm-1 


0  R,  ...  Rm 
1  m 


Calculate  k. 


m*i 


4.  Define  a  matrix  0(kj) 


1  -k 


-k  1 


5.  Produce  new  generator  matrix 

T  J 
6  -  8(k. )  G 
1  *  1 


6.  Go  to  step  2  and  repeat  until  all  ks  are  found. 

We  have  used  the  above  algorithm  to  calculate  k-values. 

Equations  (5)  and  (6)  can  be  represented  by  the  basic  lattice  diagram 
(Fig.  1).  We  shall  call  this 


Figure  1.  Basic  Lattice  Structure  (k‘+1  r  l) 


structure  the  transmission  lattice  as  will  become  obviously  justified  from  the 
following  lattice,  which  we  shall  call  the  scattering  lattice  (Fig.  2).  From 
the  scattering  lattice  it  is  clear  that  the  backward  error  and  forward  error 
can  be  viewed  as  right  going 


Figure  2.  Scattering  Lattice  Structure  (  kffl+1) 


and  left  going  signals,  respectively;  and  the  transmission  lattice  relates 
quantities  at  one  port  with  those  at  the  other  port. 

If  we  drop  the  second  subscript  on  the  error  and  represent  single 
sampling  period  delays  by  z”1,  the  transmission  lattice  looks  like  this 
(Fig.  3). 


By  cascading  transmission  lattice  we  can  produce  a  "whitening"  filter  as  shown 
in  Fig.  4.  That  is. 


Figure  4.  Whitening  Filter 
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from  Fig.  4  we  have  a  relationship  between  the  outputs  5m+1  and  em+^  and  the 
single  input  of  the  following  form: 


where  the  matrix  T  is  the  overall  transmission  matrix: 


Now,  one  can  examine  the  T22  entry  of  T  for  the  characteristic  equation  or  use 
Levinson-Durbin  to  calculate  the  AR  parameters.  Before  leaving  the 
transmission  matrix  we  consider  the  case  when  error  signals  are  delayed 
according  to  a  factor  z~®,  where  8  is  a  constant.  It  follows  from 
equations  (5)  and  (6)  that  the  lattice  structure  will  be  as  shown  in  Fig.  5. 
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_Q 

Cascading  scattering  lattices  without  the  z  p  factor,  we  generate  the 
following  filter  (Fig.  6),  which  was  used  to  predict  yn  n-0,l,...,N-l.  By 
driving  the  filter  at  the  em+^  terminal  with  a  shifted  pulse,  we 


Figure  6.  Cascaded  Scattering  Lattices 

obtain  a  response  which  estimates  the  original  discrete  time  signal,  y...  If 
we  include  the  z” ®  factors  in  each  basic  lattice,  we  get  the  diagram  shown  in 
Fig.  6  but  with  z~®  delays  with  the  eQ  and  eQ  signals  treated  appropriately. 
The  z~®  factor  is  Important  in  considerations  of  EM  waves  in  layered  media. 
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WAVES  IN  A  LAYERED  MEDIUM 

We  assume  from  an  EM  field  theory  point  of  view  that  the  source  free 
medium  of  interest  is  linear  and  stationary,  but  that  it  is  divided  into 
sections  which  are  anisotropic  and  homogeneous  bounded  by  short  sections  of 
materials  which  are  isotropic  and  nonhomogeneous.  We  assume  that  right  going 
and  left  going  planar  waves  exist  in  the  different  media  which  make  up  the 
overall  medium  and  variations  in  y  and  z  directions  are  negligible.  Figure  7 
depicts  this  situation  [5]. 


1  2  m  +  1 


Figure  7.  Layered  Medium. 

The  sections  numbered  1  to  m+1  and  shaded  are  the  isotropic  —  but  not 
necessarily  homogeneous  --  boundaries.  Between  the  boundaries  are  sections  of 
anisotropic  homogeneous  materials.  To  the  right  of  boundary  m+1  and  left  of  1 
are  Infinite  isotropic  homogeneous  regions. 

Considering  the  boundaries  first,  from  transmission  line  theory,  we  can 
write  after  a  few  manipulations  the  relationship  between  right  and  left  going 
waves  as  the  following  [4]: 
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(16) 


Wr 


i  a 

"  ut7t 

-o.(x) 


1  « 
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WR 

- 1 
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For  the  Jtth  section  a^(x)  in  operator  equation  (16)  is  given  by 


1  d\ 
at(x)  "  ST 


(17) 


where  is  a  function  of  position.  We  can  solve  equation  (17)  for  rij,  i.e., 


n.(x)  -  A  exp  [2  /  a.(A)  dX] 


(18) 


where  A  is  a  constant.  Now,  it  can  be  shown  using  forward  discretization  with 

1  A  v 

step  sizes  Ax  and  At  and  for  small  -j-  factors  that  the  scattering  lattice 
structure  for  equation  (16)  is  the  following: 


Figure  8.  Transmission  line  lattice  for  boundary  l  (l  -  l,...,m+l) 
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kj  in  the  above  Fig.  8  is  a^Sx,  where  a,  is  the  average  value  of  ct^  in  the 
boundary  £.  For  a  very  small  magnitude  value  k^  is  approximately  the 
reflection  coefficient  of  EM  wave  theory.  From  Fig.  8  we  recognize  that 
Wr  m+1  is  analogous  to  im+1  and  that  W,^  m+1  is  like  em+1.  The  only  problem 
is  the  representation  of  the  delay  in  the  basic  scattering  lattice  of  Fig.  2. 

This  problem  can  be  rectified  by  considering  waves  in  regions  between 
boundaries.  We  simply  use  the  propagation  times  to  establish  the  delays. 

Since  the  medium  is  homogenous  between  boundaries,  corresponding  reflection 
coefficients  are  zeros  and  the  scattering  lattices  are  sections  of  straight 

_Q 

lines.  Then  it  follows  that  the  z  factor  (8  *  0)  of  Fig.  5  represents  the 
time  it  takes  a  left  going  wave  to  be  emitted  by  one  boundary  and  propagate  to 
the  next  boundary.  Since  from  the  basic  lattice  structure  analogous  right 
going  "waves"  must  take  longer  than  left  going  ones  (the  z~*  factor),  the 
material  between  boundaries  must  be  anisotropic  with  propagation  speed  in  one 
direction  different  form  that  in  the  opposite  direction.  So,  not  only  do  we 
have  the  lattice  models  of  Figs.  4  and  6,  we  also  have  a  layered  medium  model 
which  will  have  the  structure  of  Fig.  9. 


RESULTS 


The  above  modeling  theory  was  applied  to  data  obtained  for  experimental 
use  of  EMPEROR  where  a  22.5  cm  monopole  was  used  as  a  test  object.  The 
response  of  the  monopole  for  250  k-values  and  for  15  non-zero  k-values 
modeling  was  obtained.  The  reduced  order  (15  non-zero  k-values)  cases  result 
from  setting  statistically  insignificant  k-values  equal  to  zeros,  and 
correspond  to  elimination  of  certain  exponentially  damped  sinusoidal 
components.  Discrete  Fourier  transforms  (DFT)  and  time  estimate  signals  are 
shown  in  Figs.  10-18.  The  only  problem  in  producing  these  results  was 
selecting  the  proper  shift  for  the  input  pulse.  This  was  done  by  trying 
different  shifts  until  the  DFT  of  the  estimate  time  signal  matched  that  of  the 
measured  frequency  domain  data. 

In  Figs.  19-21  is  shown  how  well  the  scheme  performs  for  a  more 
complicated  antenna  test  object,  the  cavity-backed  spiral  antenna.  These 
values  were  reproduced  for  250  k-values,  only.  As  in  the  monopole  case  the 
lattice  model  does  yield  agreeable  results  for  frequency  domain.  For  time 
domain,  the  results  for  CB  spiral  are  not  as  good  as  those  for  the  monopole. 


CONCLUSIONS 

It  has  been  shown  that  autocorrelation  coefficients  and  reflection 
coefficients  are  key  quantities  in  structuring  lattice  models.  It  has  been 
shown  that  frequency  domain  data  can  be  used  to  build  time  domain  lattice 
models.  It  has  been  demonstrated  that  the  method  works  for  monopole  and 
cavity-backed  spiral  antenna  data. 
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Figure  20.  Reflection  Coefficient  CB  Spiral 
{250  k-values) 


